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Core-collapse supernovae are among the most fascinating phenomena in astrophysics 
and provide a formidable challenge for theoretical investigation. They mark the spectacular 
end of the lives of massive stars and, in an explosive eruption, release as much energy as 
the sun produces during its whole life. A better understanding of the astrophysical role of 
supernovae as birth sites of neutron stars, black holes, and heavy chemical elements, and 
more reliable predictions of the observable signals from stellar death events are tightly linked 
to the solution of the long-standing puzzle how collapsing stars achieve to explode. In this 
article our current knowledge of the processes that contribute to the success of the explosion 
mechanism are concisely reviewed. After a short overview of the sequence of stages of stellar 
core-collapse events, the general properties of the progenitor-dependent neutrino emission 
will be briefly described. Applying sophisticated neutrino transport in axisymmetric (2D) 
simulations with general relativity as well as in simulations with an approximate treatment 
of relativistic effects, we could find successful neutrino-driven explosions for a growing set 
of progenitor stars. First results of three-dimensional (3D) models have been obtained, 
and magnetohydrodynamic simulations demonstrate that strong initial magnetic fields in 
the pre-coUapse core can foster the onset of neutrino-powered supernova explosions even 
in nonrotating stars. These results are discussed in the context of the present controversy 
about the value of 2D simulations for exploring the supernova mechanism in realistic 3D 
environments, and they are interpreted against the background of the current disagreement 
on the question whether the standing accretion shock instability (SASI) or neutrino-driven 
convection is the crucial agency that supports the onset of the explosion. 

Subject Index: 483, 421, 423, 415, 451, 452, 242 

§1. Supernova theory in a nutshell 

Massive stars in the range between ~8 Mq and several 10 Mq develop low- 
entropy cores, in which relativistic electrons dominate the pressure. Heavy nuclei 
yield only a small, though important, contribution to providing stabilization against 
the inward pull of gravity. The core consists of the final products of the star's nu- 
clear burning history. It is surrounded by concentric shells that, from outside inward, 
contain the successively heavier ashes of all previous burning stages (Fig. 1). 

Shell burning leads to a continuous growth of the mass of the central core until 
gravitational instability finally sets in. At this time the core resembles a hot white 
dwarf close to its maximum mass of the order of the Chandrasekhar mass. It has 
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a typical diameter of about 3000 km, a central temperature around 10^*^ K (or ap- 
proximately IMeV), a central density of several lO^gcm"^, an entropy of ~0.7 to 
~1 per nucleon (A;b is Boltzmann's constant), and a proton-to-baryon ratio of 
around 0.45. The accelerating contraction and ultimately the collapse of the de- 
generate core is initiated by the shift of nuclear statistical equilibrium (NSE). As 
the temperature in the contracting core increases, high-energy photons produce a 
growing number of a particles and free nucleons. Thermal energy is thus consumed 
to overcome the binding energy of heavy nuclei. This endothermic process, which 
partially disintegrates the iron-group material that had been assembled during the 
final burning stage of the progenitor, lowers the effective adiabatic index*) below 
the critical value (~4/3) for gravitational instability. With growing density, electron 
captures on heavy nuclei and on free protons become more and more important. 
The corresponding loss of lepton number by the production and escape of electron 
neutrinos softens the pressure increase with density even further and accelerates the 
collapse (Fig. 2, top left panel). 

Neutrinos become trapped in the 
collapsing matter at a density of roughly 
lO^^gcm^^, at which conditions their 
diffusion timescale out of the core begins 
to exceed the freefall timescale of the 
gas. Despite neutrino trapping, the col- 
lapse cannot be stopped before nuclear 
matter density (~2.7 x lO^^gcm"^) is 
reached and the equation of state stiff- 
ens because of repulsive contributions to 
the nucleon interaction potential. When 
the homologously and subsonically col- 
lapsing inner core decelerates and re- 
bounds into the surrounding, super- 
sonically inf ailing layers, sound waves 
steepen into a shock front (Fig. 2, top 
right panel). This shock front expands 
into the overlying Fe-core material, but 
is quickly damped by energy losses due 
to the dissociation of Fe-group nuclei 
into free nucleons, which extracts a ther- 
mal energy of about 8.8 MeV per nu- 
cleon from the postshock matter. Only 
1-2 ms after shock formation, the veloc- 
ities downstream of the shock have thus 
become negative. Nevertheless, the shock continues to propagate outward in mass 

*-* The effective (dynamically relevant) adiabatic index is defined as the logarithmic density 
derivative of the pressure, {d\n P/dln p)^, along a fluid element's trajectory, averaged over the 
volume of the collapsing core. It governs the transition to gravitational instability and the collapse 
dynamics. 



Onion-shell structure of pre-collapse star 




Fig. 1. Schematic onion-shell structure of a 
supernova progenitor star before core col- 
lapse. Only the main elemental con- 
stituents of the different composition shells, 
which contain the products and ashes of the 
sequence of nuclear burning stages, are in- 
dicated. Note that the radial thickness of 
the layers is not drawn to scale. The red 
square indicates the inner region zoomed 
into in Fig. 2. 
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and radius because initially the high mass accretion rate of up to several solar masses 
per second leads to the accumulation of a thick layer of dense matter behind the 
shock. Only when the mass accretion rate has decayed sufficiently and the hot, 
bloated mantle of the proton-neutron star begins to shrink in response to the lepton 
number and energy loss through neutrinos, the radial shock expansion comes to a 
halt and the shock forms a stagnant accretion shock at a radius between 100 and 
200 km (Fig. 2, middle left panel). Quasi-stationary conditions apply later on with 
only slow changes of the mass accretion rate, M, neutron star mass Mns and radius 
i?ns) and neutrino emission parameters (luminosity Li, and mean spectral energy 
(ej/)). In nonexploding spherically symmetric (i.e., one-dimensional, ID) simulations 
the shock retreats and its radius follows the contraction of the nascent neutron star 
roughly according to the relation.^) 

High mass accretion rates therefore tend to damp the shock expansion while neutrino- 
energy deposition behind the shock, which depends on the product Li, (^^^*\ can 
drive shock expansion. This issue will be elaborated on further below. 

In order to successfully launch a supernova explosion, some mechanism is neces- 
sary by which the stalled shock can be revived. Such a mechanism needs to tap the 
huge reservoir of gravitational binding energy that is released during the formation 
of a neutron star. During the infall of the stellar core the energy is first converted to 
internal energy by hydrodynamic forces (i.e., compression and the viscous dissipa- 
tion of kinetic energy in matter decelerated in the accretion shock). The degeneracy 
and thermal energy of electrons and nucleons thus stored in the proto-neutron is 
subsequently radiated away by neutrinos over a timescale of many seconds. 

Deep in the highly degenerate neutron-star interior electron neutrinos, are 
first produced by electron captures on protons. On their diffusive propagation 
towards the neutrinosphere, these electron neutrinos lose some of their energy in 
absorption-reemission processes as well as in scattering reactions with electrons and 
free neutrons and protons (Fig. 3). This effect together with the gravitational set- 
tling and compression of the outer layers of the proto-neutron star initially leads to 
rising temperatures before after some seconds cooling sets in. Since the degeneracy 
is partially lifted in the hot proto-neutron star mantle, the secondary production 
of electron antineutrinos, by positron captures on neutrons becomes possible. 
Neutrino-antineutrino pairs of all three flavors are created by thermal processes, i.e., 
nucleon-nucleon bremsstrahlung and electron-positron annihilation. Pure neutrino 
reactions (Fig. 3) also contribute to the shaping of the emitted spectra of muon and 
tau neutrinos and antineutrinos {v^, i?^, Vt, I'r):^^ which are not produced by fast 
beta reactions and thus are less tightly coupled to the stellar medium. 

Even a small fraction of the huge energy reservoir of several 10^^ ergs carried 



*' The energy transfer by neutrinos scales linearly with the neutrino luminosity and the average 
interaction cross section. The latter increases roughly with the luminosity- averaged square of the 
neutrino energy. 
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Fig. 2. Schematic representation of the evolution stages from the onset of stellar core collapse 
(top left) to the development of a supernova explosion on a scale of several 1000 kilometers. 
The displayed intermediate stages show the moment of core bounce and shock formation (top 
right), shock stagnation and onset of quasi-stationary accretion {middle left), beginning of the 
reexpansion of the shock wave ("shock revival", middle right), and acceleration of the explosion 
{bottom left). Nickel formation is indicated in the matter heated by the outgoing shock, but also 
the rising bubbles of neutrino-heated ejecta and the essentially spherically symmetric neutrino- 
driven wind {bottom right) are interesting sites for nucleosynthesis. 



Core- Collapse Supernovae 



5 



away by neutrinos is already sufficient to account for the canonical explosion energy 
of a core-collapse supernova, which ranges between some lO^'^ erg to around 10^^ erg. 
It may appear astonishing that the explosion selects an energy scale that is 2-3 
orders of magnitude lower than the reservoir of available energy. The energy scale 
of the explosion, however, is not set by the neutron star binding energy but by the 
structure of the progenitor configuration. The degenerate iron core closely resembles 
a white dwarf, embedded in the more or less dense stratification of concentric stellar 
shells containing lighter nuclear burning products. The energy scale of the explosion 
is determined by the binding energy of the progenitor layers in immediate vicinity of 
the initial mass cut between proto-neutron star and supernova ejecta*). For typical 
progenitor stars the binding energy of the silicon shell and overlying layers is of 
the order of lO^'^-lO^^ ergs, similar to the binding energy of the pre-collapse Fe core. 
Any self-regulated mechanism will tend to deposit an energy of this magnitude: Once 
matter has received this amount of energy, it will tend to become unbound and will 
expand away from what is going to become the bifurcation region between compact 
remnant and supernova ejecta. 

The delayed neutrino-driven ex- 
plosion mechanisw?) is such a self- Ncutrino Rcactions in Supernovae 

regulated process, because after absorb- 
ing energy in interactions with neutri- 
nos, the heated matter expands away gg^^ processes- 
from the neutrino-heating region. This 
naturally limits the energy input. Neu- 
trinos are produced in huge numbers in Neutrino scattering: . v + A ^ v + A 
the dense and very hot interior of the 
nascent neutron star and leave the neu- 
trinospheric region with roughly black- 
body spectra and temperatures of typi- 
cally 3-6 MeV. This is considerably hot- 
ter than the gas in the layer between 
neutrinosphere and stalled shock. Neu- 
trino heating and cooling in this region Fig. 3. Summary of important neutrino reac- 
is dominated by the absorption of elec- t^^^i^ supernova core and/or nascent 

, , . I i- i • £ neutron star. ' ' The symbol v can mean 

tron neutrmos and antmeutrmos on free ^ . , 

any type or neutrino, A represents an 
neutrons and protons and the mverse of ^^^^^^ ^^^^^^^^ ^ ^^^^^ ^^^^^^^ 

these processes, or proton (p). 

Ve + ni — )-p + e~, (1-2) 
9e+P< — >n + e'^. (1-3) 

The corresponding rates per nucleon and q~ for the energy input and loss of the 
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*-* Note that the collapse roughly conserves the total energy of the infalling shells so that the 
energetic considerations are possible for the progenitor conditions. 



6 



H.-Th. Janka et al. 



stellar medium are approximately given by:*^ 



160 .i""^ , { ^ ( MeVs-i per nucleon, (1-4) 
1052 erg/s V 100 kmy J ^ > v / 

!^ 145 I . , I MeVs~^ per nucleon. (1-5) 



^2MeV^ 

Here and T^, are the luminosity and spectral temperature of either Ve or (whose 
emission characteristics during the shock revival phase are similar), the squared 
radius r measures the geometric dilution of the flux leaving the neutrinosphere, and 
r(r) is the local gas temperature. The heating and cooling rates in Eqs. (1-4) and 
(1-5) have roughly the same magnitude. This means that both processes arc in 
tight competition. Indeed, the neutrinosphere is surrounded by a cooling layer, in 
which neutrino losses dominate. Since, however, the temperature declines roughly 
like r~^, the cooling rate falls off with and therefore much steeper than the 
heating rate, which follows basically an dependence. Consequently, there must 
be a "gain radius" Rg, outside of which neutrino heating becomes stronger than 
neutrino cooling. 

Neutrinos can deposit considerable amounts of energy in the layer between Rg 
and the shock position at radius Rs, where most of the nuclei are dissociated into 
free nucleons (the preheating in the rapidly infalling, undissociated material ahead 
of the shock is small). The optical depth for and Ue absorption in the gain layer 
between Rg and Rs can be estimated to be^^ 



rgam«0.026 l^^^^^J (^O.IMq/sJ ^200kmyl lylOOkmJ ^I.BMq^ 

(1-6) 

where M and M^s are again the rate of mass accretion by the shock and the neutron 
star mass, respectively. For typical mass accretion rates of ~0. 1-0.3 Mq s~^ and 
neutrino spectral temperatures k^T^ ~ 4-6 MeV the optical depth is 0.05-0.1, which 
means that up to about 10% of the through-going neutrinos can be captured by a 
neutron or proton. 

If the neutrino energy deposition is sufficiently strong, the stalled shock can be 
revived to potentially initiate a successful supernova explosion (Fig. 2. middle right 
panel). The threshold for runaway conditions can be coined in terms of a critical 
luminosity, Ly^d^) (for v = v^. ox v^), which depends on the mass accretion rate 
of the shock. A steady accretion shock cannot be maintained when the neutrino 
luminosity exceeds this critical vahic (for a very detailed analysis, see Ref. 10)). The 
threshold value increases with higher M, because neutrino heating has to overcome 
the higher ram pressure of the infalling material. Due to its stronger gravity also a 
larger neutron star mass Mns requires a higher neutrino luminosity for shock revival. 
The dependence of Ly^c on M and Mns can be analytically estimated on grounds of 
simple considerations. Shock expansion can become strong when neutrinos heat the 
postshock matter faster than it can be advected with the accretion flow from the 
shock downward through the gain radius into the cooling layer, where the gas loses 
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its energy by neutrino emission. This requirement is confirmed by all existing 
numerical simulations (e.g., 13)-18)). Thus requiring the neutrino- heating timescale 
to be shorter than the advection timescale through the gain layer one obtains for the 
critical luminosity:^-* 

L^,c(M) OC /3-2/5 Af2/5 ^4/5 ^ (^^.7^ 

where /3 parametrizes the ratio of postshock density to preshock density*^ This 
functional dependence nicely fits the critical curves obtained in Ref. 9), but it pre- 
dicts considerably steeper power-law dependences on Mns and M than the critical 
explosion condition of Eq. (11) in Ref. 18). It should be noted, however, that in 
the latter reference the optical depth was used as a constant free fit parameter but 
should actually be expressed as function of the fundamental quantities (Mns, -M) 
that govern the overall dynamics and structure of the postshock layer. 

Convective overturn in the neutrino-heating region has been recognized to sup- 
port the onset of the explosion^^)"^^) and to lower the critical luminosity. 
Nonradial and partially turbulent mass motions do not only stretch the dwell time 
of matter in the gain layer. Convective downdrafts carry postshock material to the 
immediate vicinity of the gain radius, where neutrino heating is strongest. More- 
over, the outward rise and expansion cooling of neutrino-heated gas in buoyant high- 
entropy bubbles reduce the energy loss by reemission of neutrinos and push the shock 
farther out. The residency time of matter in the gain layer is thiis prolonged even 
more. This combination of favorable circumstances is crucial for the development of 
runaway conditions. Nonradial oscillations and sloshing motions of the shock, which 
are associated with the growth of the standing accretion shock instability (SASI^^^), 
seem to have a similarly favorable influence on the conditions for neutrino-driven 
explosions. All these effects reduce the critical value of the neutrino lumi- 
nosity in multi-dimensional models compared to the ID case.^'^^'-^^^'^''")'^^^'^^) 

When the blast wave takes off, explosive burning leads to the production of ^^Ni 
and other radioactive species in the shock- heated silicon and/or oxygen layers (Fig. 2, 
bottom left panel). Recombination of miclcons to alpha particles and Fc-group nuclei 
in the neutrino-heated high-entropy plumes contributes to the nucleosynthctic yields, 
and the neutrino-driven baryonic outflow ("neutrino wind"), which is shed off the 
surface of the nascent neutron star by neutrino heating to fill the surroundings of 
the compact remnant after accretion has ended (Fig. 2, bottom right panel), is 
considered as an interesting site for the formation of trans-iron nuclei up to A ~ 110 
in a weak r-process and of p-rich isotopes in the neutrino-proton process, depending 
on whether the wind develops a neutron or proton excess (for a review, see Ref. 27)). 
Both the electron fraction, Y^., and the entropy in the early neutrino-heated ejecta 
and in the neutrino-driven wind have a crucial influence on the nucleosynthesis, and 
both are set by the neutrino interactions of Eqs. (1-2) and (1-3). Therefore they 
depend sensitively on the emission properties (luminosities and spectra) of the V(, 
and 9(, radiated by the forming neutron star and on the expansion dynamics of the 
ejecta (which determines the time interval of intense neutrino interactions). Also 

*^ The numerical factor in the scaling relation is found to be ~(5-6) x 10^^ erg/s for /3 ~ 10, 
M = 1 Mq/s, and Mns = 1-5 M©. 



8 



H.-Th. Janka et al. 



neutrino oscillations, in particular collective neutrino flavor transformations, can 
have an impact on the neutron-to-proton ratio that develops in the ejecta on their 
way away from the neutrino source (e.g., Ref. 28)). 

It should be noted that the explosion energy of the supernova is not determined 
at the instant when the explosion sets in; in particular it is not given by the energy 
neutrinos have transferred to the postshock matter up to this time. Instead, the 
total energy (i.e., internal plus gravitational plus the initially much smaller kinetic 
energy) of the gain layer is close to zero when the shock begins to propagate outward, 
which means that this layer is only marginally unbound.^-*' This is observed in 
all numerical models of neutrino-driven explosions (e.g., 15), 16), 29)). Instead, the 
explosion energy builds up only on a longer timescale of several 100 ms to more than 
a second. There are several sources that contribute to the final explosion energy 
(see results and discussions in Refs. 15), 29), 30)). First, neutrinos continue to heat 
"cool" gas that is freshly accreted through the shock and channelled towards the gain 
radius in convcctive downdrafts to replace there the hot matter that expands in high- 
entropy bubbles driving the shock expansion. In addition, the recombination of 
free nucleons to a particles and heavy nuclei in the ejecta releases up to ~9 MeV per 
nucleon and is a very efficient source of energy for the developing explosion (0.1 Mq 
of nucleonic matter can provide up to about 1.7 x 10^^ erg of recombination energy) 
as pointed out in Ref. 29). It should be noted that the production of free nucleons 
by nuclear photodisintegration in the infalling matter — either when the gas passes 
the accretion shock or when it is compressed and neutrino heated in the accretion 
downfiows — mostly taps the gravitational binding energy and only to a smaller ex- 
tent is fuelled by neutrino absorption. The conversion of gravitational binding energy 
to nuclear photodisintegration energy during the infall is therefore an important en- 
ergy storage whose contents are released when the matter is subsequently reejected 
and cools in the explosive outflow. Similarly, the neutrino-driven wind blown off 
the hot neutron star's surface gains power from neutrino heating as well as nuclear 
recombination. Its energy can yield a significant contribution to the energy budget 
of the supernova explosion. Scheck et al.^^^ found that ~30 70% of the blast-wave 
energy can be provided by this long-term outflow from the proto-neutron star (with 
higher relative importance for stronger explosions). In contrast, nuclear burning of 
silicon and oxygen to nickel in the shock-heated outer layers yields only a smaller 
amount of extra energy for the supernova; the production of 0.1 Mq of iron-group 
material (a typical number for normal core-collapse supernovae) releases only about 
10^° erg. During the first ~l-3s seconds of the explosion all these energy souces 
have to provide enough energy to unbind the overlying stellar layers (whose binding 
energy can range between <10^° erg and about 10^^ erg, depending on the core mass 
and compactness of the progenitor star) and, beyond that, to account for the mea- 
surable kinetic energy of the supernova. If the blast-wave energy is too low, some 
of the matter swept up by the outgoing shock will not be able to escape to infinity 
and will fall back onto the compact remnant. Because much or even most of the 
nucleosynthesized ^^Ni may thus be accreted instead of being expelled, such fallback 
supernovae are expected to be faint and hard to observe. It is therefore empirically 
unclear which fraction of stellar core collapses might belong to such types of events. 
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Predictions of Signals from Supernovae 

(magneto-)hydrodynamics of stellar plasma relativistic gravity 
(nuclear) EoS neutrino physics [ progenitor conditions 



SN explosion models 



neutrinos 



gravitational waves 





lightcurves, 
spectra 



nucleosynthesis 



explosion energies, remnant masses 



explosion asymmetries, 
pulsar kicks 



Fig. 4. Input and output from supernova explosion models. The predicted signals of stellar ex- 
plosions depend on the complex interplay of a wide variety of input physics and require the 
knowledge of the initial conditions in the progenitor stars. 



Sophisticated explosion models with all necessary theoretical ingredients are es- 
sential for a deeper understanding of observable phenomena and measurable signals 
associated with the death of massive stars and therefore also for developing a better 
definition of the role of supernovae in the astrophysical context (Fig. 4). One im- 
portant aspect of the latter is the progenitor-explosion-remnant connection, which 
includes theoretical calculations of supernova energies, nucleosynthetic yields, and 
compact object masses in dependence of the properties of the progenitor stars. This 
problem is linked to the unsolved question which progenitors leave black holes behind 
instead of neutron stars. Self-consistent numerical simulations are especially needed 
when one aims for reliable predictions of neutrino signals and gravitational waves, 
which could be detected in the case of a future galactic event and would serve as 
valuable direct probes of the dynamical processes and thermodynamic conditions in 
the supernova core. As described above, also the conditions for the explosive creation 
of heavy elements are set by the interaction of the ejecta with the intense neutrino 
fluxes from the nascent neutron star. Moreover, the production of radioactive nuclei 
like ^^Ni as well as the nucleosynthetic reprocessing of the star's composition layers 
by the outgoing shock wave depend sensitively on the explosion energy and the loca- 
tion of the mass cut that separates the compact remnant from the supernova ejecta. 
Hydrodynamic instabilities during the first second of the explosion are responsible 
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for the observed large-scale explosion asymmetries that can lead to high pulsar recoil 
velocities. By deforming the supernova shock they also seed the growth of 
secondary mixing instabilities at the composition shell interfaces after the passage 
of the explosion shock. These large-scale radial mixing processes destroy the onion- 
skin structure of the progenitor and carry heavy elements with high velocities from 
the region of their formation deep in the stellar core into the helium and hydrogen 
layers and, in turn, sweep hydrogen and helium inward in velocity as well as radial 
space. This compositional mixing and the large-scale explosion asymmetries have 
important consequences for the shape of supernova lightcurves, the characteristics 
of the supernova spectra, and the time evolution of the electromagnetic emission in 
different wave bands. 

It is clear that final answers for the wealth of questions linked to all these issues 
will ultimately require a consistent as well as consolidated solution of the explosion 
mechanism of core-collapse supernovae. Admittedly, supernova theory is not yet 
there, but it has made good progress in directions that arc promising for bringing 
us closer to the goal of our efforts. Despite the incompleteness of the present under- 
standing of the problem at the heart of exploding stars, it still seems illuminative 
to explore possible implications and thus to move forward in assembling the pieces 
of a great puzzle, in which both theoretical and observational bricks need to be in- 
terwcaved. While some adopt a pessimistic point of view and concentrate on the 
empty half of the glass, lamenting about a still imperfect match of modeling results 
with measured supernova properties and proclaiming solicitousness about what they, 
ostensibly, perceive as a defocussing and deception of the field by "myths that have 
crept into modern discourse" we prefer to look at current developments from a 
more optimistic perspective and to extract motivation from the fact that the glass 
seems to be half full at least. In this spirit a recent review article^^ has highlighted 
advances that have happened over the past ten years. But also new challenges have 
emerged from the ambitious work of researchers around the globe, whose important 
contributions to various aspects of the field have refined our picture of stellar core 
collapse and the associated physical processes. 

In the following brief, focussed overview, we shall summarize some of the latest 
results obtained by the Garching group and its collaborators. Although core-collapse 
theory is in rapid flow and the field is enthusiastically rushing towards more real- 
ism by advancing the models from two dimensions to three-dimensional space, basic 
physics that plays a role in the explosion mechanism does not depend on the consid- 
ered dimension. Progress that has been achieved in the past, although necessarily 
obtained with constrained setups like axisymmetric (2D) simulations, will neverthe- 
less provide the foundations for a better understanding of the crucial ingredients in 
a working supernova mechanism. In this sense we hope that the contents of this 
article will not grow stale too fast, despite the temporary nature of many modeling 
aspects and of conclusions that can be drawn still only from a narrowed perspective. 
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§2. Generic properties of the neutrino emission 



It is now commonly accepted that the delayed neutrino heating mechanism is 
unable to yield explosions in spherical symmetry except for progenitors stars with 
0-Ne-Mg cores instead of Fe cores. These stars at the lower end of the mass range 
for supernova progenitors reach highly degenerate conditions already after central 
carbon burning so that electron captures on Mg, Na, and Ne begin to reduce the 
effective adiabatic index and enforce the collapse. Since these progenitors possess 
an extremely steep density gradient around the core, the shock accelerates outward 
in response to the rapidly decaying mass accretion rate [cf. Eq. (1-1)] and neutrino 
energy deposition powers an outflow sufficiently strong to unbind the dilute He- and 
H-layers.^^^ The energies of such explosions are, however, low (around 10^" erg, see 



A%*VVv Shen-EoS 



^.WwV+^V"+W LS-EoS 



Fig. 15 and Ref. 38)), and only little nickel is ejected (some Mq; 
which reason the supernovae must be expected to be faint. 

In all other progenitors multi- 
dimensional flows have turned out to 
be crucial for getting explosions. If 
large-scale nonsphericity, e.g. rotational 
deformation, is absent in the progeni- 
tor core, asymmetries must grow from 
small initial perturbations by hydrody- 
namic instabilities. During core infall, 
however, the conditions are not favor- 
able for such a growth. The contracting 
flow remains essentially spherical (on g 20 

large scales) up to core bounce. More- 
over, during infall electron captures on 
heavy nuclei and free protons (the for- 
mer dominate during most of the evo- 
lution) maintain a very similar struc- 
ture of the homologously collapsing in- 
ner core, and the electron fraction con- 
verges to nearly the same central value 
{Ye^c = 0.25-0.27 for a central lepton 



Ref. 39)), for 
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0.28-0.30; Ref. 43)) 



5. Inner core masses enclosed by the 
shock-formation position in iron-core- 
coUapse simulations of a large set of (solar 
metallicity) progenitor stars.*"-' The point 
of shock formation (defined by the location 
where the postshock entropy first reaches 
3 fce per nucleon) depends on the delep- 
tonization during collapse, which hardly 
varies with the mass of the progenitor. Also 
the EoS has an influence, connected to dif- 
ferent collapse timescales and composition 
differences. We used the EoS of Ref. 41) 
in comparison to the EoS of Ref. 42) with 
compressibility modulus K — 180 MeV. 



in different progenitors despite a small 
initial spread of values in the core cen- 
ter before cohapse {Yg^i « 0.425-0.445) 
and considerable differences of the ini- 
tial central entropy {s 0.6-1.2 /cb per 

nucleon). Since the mass of the homologous core scales with the instantaneous 
Chandrasekhar mass, Mjc oc Mch oc Y^, the baryonic mass enclosed by the shock 
formation radius is also very similar over a wide mass range of progenitor stars (from 
the same modeling set; Fig. 5). The shock formation position is defined here by the 
location where the postshock entropy first reaches 3 per nucleon. The degree of 
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deleptonization during collapse, however, depends on the duration of the infall until 
neutrino trapping and on composition differences (especially the free proton frac- 
tion), both of which are connected to equation of state (EoS) properties (a detailed 
study can be found in Ref. 43)). Therefore the shock formation point is somewhat 
different for different EoSs (~O.43-M0 vs. ~O.49M0 in Fig. 5). In any case, how- 
ever, the bounce shock is launched deep inside the stellar iron core, which makes 
it impossible for the shock to overcome the energy losses by nuclear photodisinte- 
gration processes in the massive, overlying shells of iron. With the small masses of 
the homologous inner core in current models, successful supernova explosions by the 
hydrodynamical bounce-shock mechanism are ruled out. 

After a short phase of maximum strength, the bounce shock is weakened again 
by the drain of nuclear energy. A negative entropy gradient emerges in the shock 
deceleration region, which quickly decays in a short (typically 30-50 ms) phase of 
"prompt postshock convection". This phase makes a gravitational wave signal^^^ 
but has only a weak influence on the neutrino emission. Somewhat later, neutrino 
heating behind the shock grows in strength because the postshock temperature drops 
and the mean energies of the radiated neutrinos continuously rise with time (Figs. 6 
and 7). A gain layer forms and a local entropy maximum builds up just outside 
the gain radius and creates the necessary (and in most cases sufficient) condition 
for the onset of convective overturn in the accretion layer behind the stalled shock. 
First Rayleigh- Taylor fingers show up typically 80-100 ms after bounce, but it can 
take roughly another '^100 ms before the accretion flow in the postshock layer be- 
comes largely perturbed by convective downdrafts and buoyant plumes. On a similar 
timescale also the SASI grows and can manifest itself (at least in two dimensional 
(2D, i.e., axisymmetric) models) by pronounced sloshing motions of the shock sur- 
face (for a detailed numerical study of neutrino-driven convection and the SASI and 
their interaction in 2D core-collapse simulations, see Refs. 17), 26)). Besides the 
convective overturn in the neutrino-heating layer, convection also develops in the 
deleptonization region inside the neutrinospheres.^^^'^^^ Associated effects on the 
neutrino emission, for example an increase of the heavy-lepton neutrino fluxes and a 
reduction of the i>e flux compared to the nonconvective proto-neutron star, however, 
are also small in the first ~ 100 ms after core bounce. 

During this early postbounce period of up to about 100 ms, spherically symmetric 
models therefore yield reliable information about the neutrino emission properties. 
We display corresponding information on the radiated luminosities and mean spectral 
energies (for an observer at rest at a large distance from the source) for all neutrino 
species in Figs. 6 and 7. (Muon and tau neutrinos and antineutrinos have very 
similar opacities and are lumped together to one kind of heavy-lepton neutrino, z^^:-) 
In the first case the relatively stiff nuclear (EoS) of Ref. 41) was used, in the second 
case the considerably softer one*^ of Ref. 42). These images should be compared 

*^ Somewhat unconventionally, we use the terms "stiff" and "soft" for the EoS to classify the 
radius evolution of the forming neutron star with the corresponding EoS, "soft" meaning that the 
nascent remnant contracts faster to its final radius. Note that a "soft" EoS in our sense neither 
means that the maximum neutron star mass is low nor that the final neutron star radius is partic- 
ularly small. Physically, the maximum mass of the neutron star is determined by the supernuclear 
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Fig. 6. Time evolution of luminosities (/e/f) and mean energies (defined as ratio of energy flux to 
number flux; right) for (top), Ve {middle), and one kind of heavy-lepton neutrino (all treated 
equally; bottom) as measured in the lab frame at infinity. Time is normalized to the moment 
of core bounce. The results are plotted for (nonexploding) ID simulations of a set of solar- 
metallicity progenitors*''' and an older ISM© progenitor (sl5s7b2; Ref. 46)), applying the EoS 
of Ref. 41). 



properties of cold neutron star matter, whereas the radius evolution of the nascent compact object 
strongly depends on how the conditions in matter at subnuclear densities depend on changes of the 
temperature and of the neutron-to-proton ratio (both are evolving in response to neutrino emission). 




-0.1 0.0 0.1 0.2 0.3 0.4 -0.1 0.0 0.1 0.2 0.3 0.4 
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Fig. 7. Same as Fig. 6, but for (nonexploding) ID simulations with the EoS of Ref. 42) using a 
compressibihty modulus of K — 220 MeV. This EoS is softer than the one used in Fig. 6. 



with Fig. 3 in Ref. 47). Although our set of core-collapse calculations and the one of 
Ref. 47) employ different series of progenitor models, basic features of the neutrino 
signal are generic and, even more, general trends should be comparable even beyond 
100 ms after bounce, despite the disregard of multi-dimensional effects. 

It should be noted that our calculations were performed with the Prometheus- 
Vertex code with general relativistic corrections,^^ which performed very well in 
tests against fully relativistic treatments. ^^■^^•''^^^ The energy-dependent three-flavor 
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neutrino transport module is based on a two-moment closure scheme with a variable 
Eddington factor derived from a model-Boltzmann equation. The transport includes 
velocity-dependent terms, energy-bin coupling, and the full set of neutrino opacities^^ 
(see Fig. 3). 

Although in Ref. 47) a variety of simplifications were applied (e.g., a subset 
of neutrino processes was used and velocity-dependent effects were ignored in the 
transport), the overall evolution of the neutrino emission properties of Ve and Ve is 
reproduced relatively well. This outcome is connected to the facts that on the one 
hand electron-flavor neutrinos mainly interact by the charged-current absorption and 
emission processes and on the other hand the neutrinos escape from regions near the 
proto-neutron star surface, where the matter is essentially at rest. Features that 
depend on the motion of the stellar medium, however, are absent in the results of 
Ref. 47). One example is the characteristic local minimum of the U(, luminosity 
shortly after core bounce and before the luminosity rises steeply to its prominent 
maximum. In this feature the luminosity drops from ~1.15 x 10^^ ergs~^ to ~0.7 x 
lO^^ergs"^ (which can easily be seen only on a zoom to its location). This dip is 
formed about 1 ms after bounce and approximately another millisecond before shock 
breakout from the neutrinosphere (where the optical depth decreases to about unity). 
At this time the shock is still deep in the optically thick region and the neutrinos 
produced behind the shock are still trapped. Most of the release therefore comes 
from the compressed, unshocked and semitransparent layer around 100 km. When 
the matter in this layer falls inward with increasingly higher speed, the emission gets 
more and more redshiftcd due to the Doppler effect and, moreover, an increasing 
fraction of the neutrinos is advected inward with the gas flow. These effects lead to 
the transient luminosity reduction before the rise to the shock-breakout burst sets 
in. 

The U(, signal through core collapse, bounce, and shock breakout is extremely 
similar for all progenitors and also for different EoSs.^'^^ Even during the first 
~ 100 ms after bounce the progenitor stars have little influence on the mean en- 
ergies of the radiated neutrinos of all kinds, although the luminosities already begin 
to exhibit a spread that reflects the different mass infall rates in the cores of different 
stars. More compact progenitor cores*-* lead to higher mass accretion rates of the 
nascent neutron star, which in turn produces higher neutrino luminosities supplied 
by the emission of neutrinos from the hot, lepton-rich, freshly accreted material. 
Higher stellar core compactness systematically correlates with higher luminosities of 
all neutrino species.^^^ 

At times t > 100 ms post bounce, consequences of multi-dimensional effects on 
the neutrino emission cannot be ignored any longer.^^^'^^^'^^) The radiated neutrino 
signal is then affected by convection inside the neutron star, which raises the heavy- 
lepton neutrino fluxes, reduces the Ug emission, and decreases the mean energies 
of the escaping neutrinos. ^^^'^^^ The neutrino emission is also modified and mod- 



*' Compactness is measured by the ratio of mass to corresponding radius enclosing this mass. 
Higher compactness means that more mass falls inward in a certain time. For a definition and 
discussion of the compactness in the context of black hole formation, see Ref. 51). 
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ulated by nonspherical accretion downflows, which are associated with anisotropic, 
transient, short "bursts" of accretion luminosity. ^^^"^^^ Nevertheless, some general 
properties can be discussed that are shared by ID as well as angle-averaged 2D and 
3D results. 

At i > 100 ms post bounce also the mean energies of the radiated neutrinos begin 
to exhibit a spread that reflects the core compactness. Higher core compactness and 
higher mass accretion rates do not only produce higher neutrino luminosities but 
also larger mean energies of the emitted neutrinos. While before t ~ 100 ms p.b. the 
nascent neutron stars in all progenitors are still relatively similar, the different mass 
accretion rates in stars with different core compactnesses cause the proto-neutron 
star masses to become increasingly different. For this reason the mean energies 
develop a growing spread with time: More massive proto-neutron stars possess hotter 
neutrinospheres and the mean neutrino energies increase more steeply than in the 
case of less massive neutron stars in less compact progenitors. The collapsing cores 
of such less compact stars feed the growth of the neutron star mass with rapidly 
decreasing rates at late postbounce times. Interestingly, in all cases (and for both 
stiff and soft neutron star EoS) there is a point during the late accretion phase when 
the mean energy of Ug crosses and subsequently exceeds the mean energy of the 
heavy- lepton neutrinos. The physical reason for this effect was discussed in Ref. 52). 
After the explosion has taken off and the accretion has ended, the mean energies 
adopt again the well known hierarchy, where {cu^} < {epj < (e^^). 

It should be noticed that in Ref. 47) — presumably because of the approximations 
in the neutrino transport used there — considerably higher mean neutrino energies 
are predicted for the later postbounce evolution (compare Fig. 3 in Ref. 47) with 
Figs. 6 and 7 in our paper*)). In particular, also the differences between the core- 
collapse models with soft and stiffer EoS tended to be larger in Ref. 47). The 
authors of the latter paper proposed to use the cumulative neutrino energy from 
the postbounce accretion phase, which correlates with the mass accretion rate of 
the progenitor, to obtain information of the core structure of the collapsing star 
in the case of a neutrino measurement from a future galactic supernova. In order 
to break the degeneracy between nuclear EoS and progenitor-core compactness (a 
higher neutrino luminosity and number flux, for example, can be caused by a larger 
mass accretion rate or by a softer EoS), they suggested to combine the measurement 
of the cumulative energy (or cumulative event number) with the information about 
the mean energy of the neutrinos. In view of the weak sensitivity of the mean energy 
to the EoS (cf. Figs. 6 and 7) such a discrimination will be quite an ambitious 
undertaking. 



*' Note that the comparison between the different transport solvers in Fig. 2 of Ref. 47) was done 
with a reduced set of neutrino opacities instead of the more sophisticated and larger set of neutrino 
reactions usually applied in the simulations of the Garching group. Moreover, the approximate 
description of general relativistic gravity used in the Prometheus- Vertex calculations of Ref. 48) 
(referred to in Ref. 47)) was subsequently improved by Marek et al.,"'*-' whose preferred description 
was shown to lead to much better agreement with fully relativistic results.®^ 
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Fig. 8. Pre-coUapse density profiles of pro- 
genitor stars used for the 2D explosion 
models of Figs. 10-15. The circles mark 
the boundaries of the Fe-core and the 
Si/Si+O interface in the solar-metallicity 
11.2 M© (sll.2), 25 Mq (s25.0), and 27 Mq 
(s27.0) models,*"^ in the 15 Mq (sl5s7b2) 
model, and in the ultra-metal poor 
(10"" solar) 8.1 Mq progenitor,^®' whereas 
they indicate the Fe/Si and Si/O+N tran- 
sitions in the zero-metallicity 9.6 Mq (z9.6) 
star.^®' In the 8.8 Mq progenitor with O- 
Ne-Mg core^®' the circles denote the com- 
position changes from NSE to Ne-I-O, then 
to Ne-I-Mg-I-O, to C-l-O, and finally to He. 
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time after bounce [s] 

Fig. 9. Postbounce evolution of the maximum, 
minimum and average shock radius in a 
2D simulation of the 27 Mq progenitor 
of Fig. 8 with the Prometheus- Vertex 
code, which employs an approximate treat- 
ment of general relativistic effects in the 
gravitational potential and neutrino trans- 
port. As in the fully relativistic 2D simula- 
tion shown in Fig. 10, bottom right panel, 
a neutrino-driven explosion takes place, al- 
though the runaway shock expansion sets 
in somewhat later and after less vigorous 
bipolar shock oscillations than in the rela- 
tivistic calculation. 



§3. Relativistic explosion models in two dimensions 

Successful neutrino-driven explosions cannot be obtained in ID with the neu- 
trino heating associated with the neutrino emission properties of present models as 
discussed in Sect. 2. The Garching group, however, has recently found explosions 
for a growing set of progenitor stars with different masses (8.1, 8.8, 9.6, 11.2, 15, 
and 27 Mq) and metallicities (Fig. 8) in 2D simulations with the general relativistic 
CoCoNuT- Vertex hydrodynamics and neutrino transport code.^^'"^^)'"^'^)'^^) These 
results can be considered as a basic confirmation of the self-consistent 2D explosion 
models of 11.2 and ISM© stars computed with the Prometheus- Vertex code 
by Buras et al.^^-* and Marek & Janka.^^^ Also the 27 Mq explosion has recently 
also been reproduced (Fig. 9) with the Prometheus- Vertex tool, which employs 
a relativistic approximation of the gravity potential^^^ and relativistic corrections in 
the transport module. Despite differences in details, e.g. in the postbounce dy- 
namics and explosion times of 2D models with general relativistic and approximate 
relativistic treatment, there seems to be overall agreement and compatibility of the 
results obtained with these two different numerical schemes*^. 



*' The Prometheus- Vertex and CoCoNuT- Vertex codes do not only differ in the treatment 
of general relativity, they also employ largely different hydrodynamics solvers, cf. Refs. 4) and 6). 
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Fig. 10. Time evolution after core bounce of the profiles of entropy per nucleon in north-polar and 
south-polar directions for 2D (general relativistic) explosion simulations^®'' ^'^■''^^^ of the 8.8 
ONeMg-core progenitor and the 8.1 M©, 9.6 Mq, 11.2 Mq, ISM© and 27 Mq Fe-core stars with 
different metallicities {from top left to bottom right) shown in Fig. 8. Dark blue and black 
represent low-entropy unshocked matter while green, yellow, and red indicate hot, high-entropy 
matter. Bipolar shock oscillations suggest SASl activity. 



Entropy profiles in the north and south polar directions as functions of post- 
bounce time are displayed for the relativistic explosion models in Fig. 10. Corre- 
sponding Ye and entropy distributions for representative instants near the end of 
the simulated postbounce periods are given in Fig. Of the set of progenitors 

investigated so far, only the 25 Mq did not show any tendency of an explosion and 
looks very disfavorable for a final success when the computational run was stopped 
at 450 ms after bounce (Fig. 12). This can be seen in Fig. 14, which shows a very 

*-* Note that the explosion simulation for the 8.8 Mq progenitor with O-Ne-Mg core was per- 
formed with the relativistic approximations of the Prometheus- Vertex scheme. 
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Fig. 11. Cross-sectional distributions of electron fraction (left half panels) and entropy per nucleon 
(right half panels) for the simulations shown in Fig. 10 {from top left to bottom right for 8.8, 
8.1, 9.6, 11.2, 15, and 27 M© progenitors) at postbounce times near the end of the simulations: 
tph = 365 ms, 330 ms, 318 ms, 920 ms, 775 ms, and 790 ms, respectively. The supernova shock is 
marked by a thin white line. High-entropy bubbles of neutrino-heated, expanding matter drive 
aspherical shock expansion in most cases. 
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Fig. 12. Same as Fig. 10, but for the non- 
exploding solar-metallicity 25 Mq star of 
Fig. 8. Alternating north-polar and south- 
polar shock excursions indicate violent 
SASl sloshing motions. 



0.55 



0.46 



0.37 



0.28 



0.19 



0.10 






X [km] 



24 
20 
16 
12 
8 
I 4 


50 100 

5 [A:/nucleon] 



Fig. 13. Same as Fig. 11, but for the nonex- 
ploding solar-metallicity 25 Mq model of 
Fig. 12 at 309 ms post bounce. The post- 
shock layer is very narrow and the rapid 
infall velocities of the fluid between shock 
and gain radius suppress the growth of 
neutrino-driven convection. The entropy 
structures in the postshock layer are mainly 
linked to Kelvin-Helmholtz instability asso- 
ciated with strong shear effects at the inter- 
face of differentially moving regions in the 
sloshing flow. 



small maximum radius of the shock in the 25 Mq case, while in all other relativistic 
2D models the shock expands with high velocities of at least 6000-7000 km s~^, in 
most cases much faster. The explosions develop on different timescales, depending 
on the different density profiles and locations of the composition shell interfaces, 
which determine the decay of the expansion-damping mass accretion rate. (Remem- 
ber that according to Eq. (1-7) the critical luminosity is larger for higher M and 
bigger neutron star mass, which also grows faster for high M.) The decrease of the 
mass accretion rate with time competes with the progenitor-dependent evolution of 
the neutrino emission properties (cf. Figs. 6 and 7), which decide about the strength 
of the neutrino energy deposition. If the critical condition for an explosion in 2D is 
met at some point, the model makes the transition to a runaway expansion of the 
supernova shock. In the 25 Mq star with its highest compactness (Fig. 8) the large 
mass infall rate (and quickly growing neutron star mass) prevent an explosion until 
the end of our computation. 

In Fig. 15 the "diagnostic energies" are given for the successful cases. The diag- 
nostic energy at each time is defined as the total energy (internal plus gravitational 
plus kinetic) of all postshock material with positive radial velocity and and positive 
total specific energy. This is not yet the explosion energy of the supernova, because 
the positive energies may still increase by ongoing neutrino heating of accreted and 
then reejected gas, by additional energy associated with the neutrino-driven wind 
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Fig. 14. Maximum shock radii versus time for 
the set of 2D explosion models of Figs. 10 
and 11 and for the nonexploding 2D model 
of Figs. 12 and 13. 
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Fig. 15. Evolution of the "diagnostic" explo- 
sion energy (i.e., the energy of outward 
expanding postshock matter with positive 
specific energy) versus time for the set of 
models of Fig. 14. 



material, by the recombination of free nucleons and a particles to heavy nuclei, and 
by nuclear burning in the shock-heated layers (see Sect. 1). Moreover the binding 
energy of the outer stellar layers (still ahead of the shock) will contribute on the 
negative side and has to be overcome to make the star unbound for an explosion 
with excess kinetic energy at infinity. While this binding energy is negligible for the 
8.8 Mq progenitor (whose 0-Ne-Mg core is surrounded by an extremely dilute and 
loosely bound H-shell) and small for the lower-mass iron-core stars (i.e., less than 
~3 X 10^^ erg in the 8.1 case and even less for the 9.6 Mq star), the binding 
energy of the stellar mantle and envelope can be appreciable for more massive pro- 
genitors. In the 11.2 M0 and 15 Mq models, for example, the binding energies of 
the preshock shells are 7.5 x 10^^ erg and 2.6 x 10^" erg, respectively, whereas the 
still available recombination energies of the postshock matter are 2 x 10^^ erg and 
1-2 X 10^*^ erg. In the case of the 15 Mq and 27 Mq models the diagnostic ener- 
gies are still steeply increasing at the end of the simulations and the terminal values 
cannot be guessed. In contrast, the diagnostic energies of the other models begin to 
saturate. The explosions will therefore become relatively weak, and a lot of fallback 
has to be expected. 

§4. Magnetically supported explosions in nonrotating progenitors 

Magnetic fields are known to be strongly amplified during the secular postbounce 
evolution of rapidly differentially rotating, collapsing stellar cores. Besides compres- 
sion upon infall, the amplification is mainly achieved by field winding in shear layers 
and the magnetorotational instability (MRI), see e.g., Refs. 57), 58). One may ask, 
however, how strong the fields have to be, and how strong the corresponding initial 
fields, if magnetic effects are to have an impact on the development of the explosion 
in nonrotating (or very slowly rotating) stellar cores. Slow core rotation is predicted 
by evolution models of magnetized, massive stars^^^ because of angular momentum 
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loss mainly during the red-giant phase. With predicted pre-collapse spin periods of 
more than hundred seconds in the iron core, the postshock layer attains rotation pe- 
riods of hundreds of milliseconds to many seconds and possesses only tiny amounts 
of free energy of rotation (which is the reservoir to be tapped for field amplification). 
In this case neither wrapping nor the MRI are efficient mechanisms to enhance the 
field strength. 

Obergaulinger & Janka^'^^ have re- 
cently investigated the effects of mag- 
netic fields on the shock evolution in 
nonrotating stellar cores*\ In this case 
the amplification of the fields in the gas 
flow relies on two processes: 

• Compression of the fields against 
the magnetic pressure increases 
the magnetic energy density emag; 
magnetic energy is created by this 
process at a rate Scmp = — emagV-v. 

• Stretching and folding of the field 
lines, in which case the energy den- 
sity of the field can be amplified at 
a rate of Sstr = B'^B^V jVi. This 
process is the main ingredient of a 
turbulent dynamo. 

The amplification by compression takes 
place in the converging flow of the un- 
shocked, infalling layers, in the flow 
passing the shock, and in convec- 
tive downdrafts in the postshock re- 
gion. Turbulent mass motions due 
to neutrino-driven convection and the 
SASI^^^'^^) lead to further ampliflcation 
by dynamo action in the region between 
stalled shock and protoneutron star. 

Figure 17 shows the postbounce 
evolution of four simulations of a 15 Mq 




Fig. 16. Field lines and field strength (color- 
coded) with central value of Be = 10^'' G 
assumed for the initial (pre-collapse) field 



configuration of a 15 Mq progenitor. 
Initial configurations with other central 
field values used in our simulations are ob- 
tained by simple scaling. In our nonrotat- 
ing stellar cores we have adopted a purely 
poloidal field. This is clearly artificial and 
besides the experimental purpose the as- 
sumed absence of stellar rotation to cre- 
ate toroidal field components might serve 
as justification. A high-resolution version 
of this figure is available upon request. 



progenitor, comparing the nonmagnetic 

case with three cases for different initial magnetic fleld strengths, corresponding to 
central fleld values of 10^^ G, 10^^'^ G, and 10^^ G. The assumed pre-collapse fleld 
conflguration in the stellar core is displayed in Fig. 16. The panels of Fig. 17 visualize 
the time dependence of the radial proflles of angularly averaged entropy and magnetic 
fleld values (or in the nonmagnetic case). Spatial distributions for representative 



*•* The results discussed here are based on simulations with higher resolution and a two- 
dimensional, multi-energy-group, two-moment closure scheme for neutrino transport, which was 
upgraded relative to the models in Ref. 60). The improvements will be entered in a revised version 
of the latter paper. 
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Fig. 17. Postbounce evolution of 2D simulations for a 15 Mq progenitor'*"-' with different assump- 
tions of the initial magnetic field strength for central values of _Bc = {top left), 10^* G {top 
right), 10^*'^ G {bottom left), and 10^^ G {bottom nght). In the case of the nonmagnetized 
model, radial profiles of the angular averages of entropy per baryon (upper half panel) and 
electron fraction, Ye, are color-coded, whereas for the magnetized models Fg is replaced by the 
logarithm of the angularly averaged magnetic field strength. Note that the ordinate with the 
radius has been mirror imaged in the lower half-panels. The thin grey lines indicate trajectories 
for different fixed values of the enclosed mass spaced in steps of 0.025 Mq (dashed) and 0.1 Mq 
(solid). The two black lines denote the maximum radius of the shock front and the outer edge 
of the convection zone inside of the proto-neutron star, respectively. The white dashed lines 
mark locations of chosen constant values of the angle- averaged mass density for p = lO^"*, 10^^, 
lO^^gcm"^, etc. from the center outwards. The amplification of the initial core field during 
collapse, shock compression, and turbulent folding fosters shock expansion and the initiation of 
neutrino-powered explosions on a timescale that decreases with the strength of the magnetic 
field (compare left lower and right lower panels). 
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Fig. 18. Representative snapshots for the simulations displayed in Fig. 17. For the nonmagnetized 
model the radial velocity (left half panel) and the entropy per nucleon (right half panel) are 
shown; for the models with magnetic fields the logarithm of the magnetic field strength is color- 
coded instead of the radial velocity. Arrows indicate the flow field. The chosen postbounce 
times are 400ms, 512ms, 574ms, and 412ms (from top left to bottom right), respectively. 
Neutrino-heated, high-entropy gas rises in buoyant bubbles and pushes the shock expansion. 
The magnetic field structure in the postshock layer traces the convective pattern of bubbles and 
downflows and tends to be stronger in the compression regions of the downdrafts. 



times are presented in Fig. 18. 

In the lower two panels of Fig. 17 the shock evolution is clearly affected by 
the presence of the magnetic fields. The shock radius is larger compared to the 
nonmagnetic model (upper left panel) and ultimately the shock begins to expand, 
followed by the gas in the postshock region. While for an initial central field strength 
expansion becomes strong about 200 ms after bounce (bottom 
right), a powerful shock expansion sets in only 500 ms after bounce for initial central 
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Fig. 19. Postbounce evolution of the 2D simulations with initial central magnetic fields of 10^^'^ G 
[left) and 10^^ G [right). The upper half-panels show the color-coded logarithm of the ratio 
of magnetic pressure to gas pressure, the lower half-panels display the logarithm of the ratio 
of magnetic energy density to kinetic energy density. The shock expansion is supported by 
significant magnetic pressure (between several per cent and ~10% of the gas pressure) in the 
postshock layer. The field there reaches equipartition strength with the kinetic energy of the 
flow. 



field strength of 10^^'^ G (bottom left). For even weaker initial field {Be = 10^^ G; 
upper right panel) one can see first hints for a transition to shock expansion after 
nearly 600 ms. The onset of shock acceleration is associated with the development 
of positive total energy in the expanding postshock material. This growth of the 
energy signals the beginning of an explosive runaway. 

The dynamical relevance of the magnetic field can be expressed by the ratio 
of magnetic field pressure to the gas pressure, whose radial profiles are plotted as 
function of time for the two exploding models with highest fields in the upper half- 
panels of Fig. 19. Postshock fields of about 10^^ G (Fig. 17) correspond to a magnetic 
pressure contribution to the gas pressure of several per cent up to about 10%. This 
obviously is sufficient to play a dynamically important role and to support the shock 
expansion. As a consequence, the inflation of large bubbles of neutrino-heated, high- 
entropy gas is enabled and the flow in the accretion region between shock and gain 
radius becomes much more ordered, with a dominance of low-order spherical har- 
monics modes in the convective pattern (Fig. 18). As the shock propagates outward, 
the push of the rising bubbles of heated matter becomes the driving force for fur- 
ther expansion, because the relative importance of the magnetic pressure to the gas 
pressure behind the outward rushing shock quickly decreases. The lower half-panels 
of Fig. 19 demonstrate that in the turbulent postshock layer prior to the onset of 
the explosion the average magnetic field grows to a value near equipartition with the 
amplification-supporting average density of kinetic energy. 

Despite the constraint of our models to 2D and the unavoidable limitations of 
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the resolution, we think that one can conclude that the magnetic field amplification 
in collapsing stellar cores, even in the absence of any significant amount of initial 
core rotation, can lead to dynamical effects on the evolution of the supernova shock. 
For sufficiently strong initial fields magnetic pressure in the postshock layer can 
support the onset of neutrino-powered explosions on relevant timescales (of order 
0.5-1 s after core bounce). Further simulations, in particular also 3D models with 
more realistic initial conditions for the field strength and geometry based on stellar 
evolution predictions, are needed to clarify whether this could play a relevant role 
as ingredient in the neutrino-driven explosion mechanism. 

§5. Convection or SASI as trigger of shock expansion? 

A detailed discussion and analysis of the postbounce dynamics leading to the 
onset of explosions in the 2D models of Sect. 3 (including a description of the associ- 
ated gravitational- wave signals) is provided in Rcfs. 15)-17),44). Some of the plots 
of the entropy evolution in Figs. 10 and 12 as well as corresponding movies show 
violent bipolar sloshing motions of the shock and of the matter in the postshock 
layer. 

These bipolar shock oscillations are found to exhibit growing amplitudes before 
neutrino-powered explosions set in. Such a behavior is visible in the 11.2 M0 and 
15 Mq cases, ^^^'^^^ but the alternating up and down motions of the shock are espe- 
cially prominent in the 27 Mq model, where the extremely high infall velocities 
in the postshock flow suppress the growth of neutrino-driven convection in the first 
place, unless large initial entropy perturbations drive non-linear buoyancy from the 
beginning (for a discussion of the necessary conditions, see Refs. 13), 17), 26), 63)). 
In this case the SASI-typical oscillatory increase of the amplitude of low-order shock 
deformation modes appears very clearly. In contrast, the dynamics of the lower- 
mass progenitors does not exhibit this characteristic behavior or, at least, it is not 
present there in such a pure form (Fig. 10). Buras et al.^^^ found shock oscilla- 
tions with large dipole {£ = 1) amplitudes to be crucial for the explosion of the 
11.2 Mq case because they obtained successful shock revival in 180-degree (pole-to- 
pole) simulations, whereas simulations with a 90-degree wedge around the equator 
did not explode.*). The violent shock motions and the associated shock expansion 
allow matter accreted through the shock to stay in the neutrino-heating region for 
a longer time. The increased abidance timescale in the gain layer enhances the ef- 
ficiency of neutrino-energy transfer and thus supports the development of runaway 
conditions. 

The shock-oscillation phenomenon exhibits strong similarity to the SASI dy- 
namics that can be observed in 2D as well as 3D simulations of adiabatic accretion 
flows using an ideal-gas EoS and simple representations of the effects of neutrino 
cooling as regulator of (quasi-)stationary postshock conditions.^^^'^^^'^^) An exper- 
imental shallow-water analog that considers a hydraulic jump in a converging two- 



*' Successful shock revial, however, could also be obtained for the 11.2 Mq star with a 90-degree 
pole-to-equator grid. 
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dimensional water flow has also been developed®^) to demonstrate the basic aspects 
of the instability that preferentially leads to low-order (dipolar, quadrupolar) shock 
deformation and violent, nonradial shock expansion and contraction. While the axis 
of these motions naturally coincides with the artificial symmetry (polar-grid) axis 
in 2D-axisymmetric simulations, there is no such predetermining constraint in the 
experimental setup. This confirms that the grid geometry is not the crucial aspect 
that allows the SASI sloshing phenomenon to occur. Moreover, the growth behav- 
ior of the SASI from the linear regime to the nonlinear stage in competition with 
neutrino-driven convective instability was investigated for supernova core conditions 
in 2D by Foglizzo et al. and Scheck et al. in Refs. 26),63)*). In the linear regime 
both instabilities can be discriminated by their different growth behavior (oscilla- 
tory vs. nonoscillatory) and by the absence or presence of buoyant flow structures in 
the postshock shell. Stronger neutrino heating or slower advection velocities in the 
accretion flow were recognized to be favorable for the growth of neutrino-driven con- 
vection, whereas for insufficient heating and in fast accretion flows the SASI modes 
were seen to grow predominantly. Convective activity can also win when its growth 
is accelerated by strong buoyancy forces acting on large initial entropy perturbations. 
In the fully nonlinear regime, however, a clear separation of both instabilities could 
not be achieved and is very difficult, because SASI shock motions trigger secondary 
convection,^^) and inversely, it is also conceivable that neutrino-driven convection 
feeds back into enhanced SASI activity. Fundamental aspects of these hydrodynam- 
ical results seem to be in agreement with global linear stability analysis of stalled 
shocks in accretion flows (using a microphysical EoS and taking neutrino reactions 
into account) by Yamasaki & Yamada,^^-* who found the dominant growth of oscil- 
latory or nonoscillatory radial and nonradial modes, depending on the size of the 
imposed neutrino luminosities. 

For all these reasons Marek & Janka^^^ (see also Ref. 13)) interpreted the violent 
shock-sloshing motions with growing amplitudes, which preceded the explosions they 
found in 11.2 Mq and ISM© simulations, as a consequence of SASI activity. They 
were, however, fully aware that this was a speculative interpretation, a working 
hypothesis, which was by no means based on unambiguous facts and well-founded 
theoretical arguments. A solid theoretical foundation was (and still is) lacking for the 
nonlinear stage of the observed phenomenon. Nevertheless, it is clear that the energy 
input for the explosion has to be delivered by neutrino heating (whose behavior was 
quantitatively analysed and energetically evaluated by Marek & Janka) , and that the 
shock oscillation motions are only helpful for improving the heating conditions, but 
by themselves are not the driving agency of the explosion. In this context it should 
be kept in mind that the growth of the SASI depends on an advective-acoustic 
cycle. Once shock expansion initiates mass expansion behind the shock as well, 
the cycle tends to lose its support because the inward flow of matter is decelerated, 
diminished, or even quenched. It therefore seems unlikely that SASI is the dominant 



*' Numerous analytical and numerical studies for the linear growth regime of the SASI in su- 
pernova core like environments been performed. For probably the most comprehensive discussion 
including relevant references, see Ref. 67). 
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mode of nonradial mass motions up to the stage where the explosion takes off. It may, 
however, play a supportive role for establishing favorable neutrino-heating conditions 
and for getting close to the final runaway situation. 




Fig. 20. Volume- rendering images of two instants (180 ms and 292 ms post bounce) of a first ex- 
plorative 3D simulation of an 11.2 Af0 progenitor (model sll.2 of Ref. 40)) with full energy- 
dependent ray-by-ray three-flavor neutrino transport and the sophisticated set of neutrino opac- 
ities as applied also in the ID and 2D models of the Garching group shown in Figs. 5-7 and 
9-15. Neutrino- heated bubbles dominate the asymmetries in the postshock layer. The absence 
of large-amplitude bipolar SASI shock motions in this simulation challenges speculations in the 
literature''^'' that local neutrino heating computed with the ray-by-ray approximation might 
unphysically correlate with low-order mode shock and matter motions and could be pumping 
them unphysically. 

Burrows^^-* classified the SASI as a "myth that has crept into modern discourse" . 
(In fact, it has sneaked into his own works, too, cf., e.g., Refs. 53), 69).) In his view 
the vigorous dipolar shock oscillations seen in 2D simulations with full neutrino trans- 
port have been incorrectly associated with the SASI sloshing motions observed in 
simplified, neutrino- free studies. He argued that due to the inverse energy cascade of 
2D turbulence, which transports energy from small to large scales, neutrino-heated, 
buoyant bubble convection with dominant low-^ modes has been confused with the 
SASI, whose growth rates are highest on the largest angular scales (i.e., for the spher- 
ical harmonics modes of lowest order i > 1). Such a possibility cannot be rejected 
categorically, it is indeed a viable possibility. Just as neutrino energy input might 
play a role as driving force of enhanced SASI activity, neutrino-driven buoyancy 
could also play a destructive role by destroying the flow coherence needed for SASI 
ampliflcation. 

But a proof of a negative feedback of neutrino-driven convection on the SASI 
has so far not been given, and a convincing demonstration of an exclusive dominance 
of neutrino-triggered buoyancy has also not been provided. Burrows et al.,^''^ for ex- 
ample, argued that in their two-dimensional numerical toy model stronger neutrino 
heating and thus enhanced neutrino-driven convection was the cause of the growth 
of the dipolar sloshing amplitude that was crucial for getting neutrino-driven explo- 
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sions. This result, however, may be tightly connected to the simplified setup chosen 
for their simulations, especially with respect to the treatment of the neutrino source 
terms, by which energy losses from the cooling layer near the neutron star surface 
were massively underestimated. Moreover, neutrino emission from the neutron star 
interior was ignored and therefore the contraction of the neutron star was also under- 
estimated. Investigations based on such very particular conditions without varying 
the progenitor structure, neutron star mass, gravity, and contaction behavior, as well 
as the efficiency of neutrino cooling, are endangered to lead to conclusions that hold 
for this special case but do not necessarily possess wider validity. 

In consequence of the mentioned deficiencies of the neutrino and neutron star 
treatment, the accretion velocities are reduced, which implies more favorable condi- 
tions for the growth of neutrino-driven buoyancy while at the same time a possible 
importance of the SASI is suppressed, because the SASI is amplified more rapidly 
for shorter advection timescales in faster accretion flows (see, e.g., Ref. 26)). It is 
therefore not astonishing that Burrows et al.^'^) observed only small SASI amplitudes 
in the absence of neutrino heating and the associated buoyant convection. It is also 
little astonishing that Murphy et al.^^^ (also Ref. 72)) diagnosed the conditions in the 
same simplified modeling setup as fully consistent with neutrino-driven convection 
and turbulence in 2D as well as 3D simulations. The 27 Mq calculation discussed 
above and in Sect. 3, which includes a sophisticated and self-consistent treatment 
of all relevant physics aspects, demonstrates that there are situations in collapsing 
stellar cores at least under the constraint of axisymmetry and most obvious in 
the 25 Mq and 27 Mq cases — where the SASI is the most quickly growing insta- 
bility and initially drives the shock expansion. It thus establishes, of course, also 
better conditions for enhanced neutrino heating. This may at some point change 
the situation from SASI-dominated to buoyancy-dominated, because it is clear that 
increasingly stronger neutrino heating will foster stronger buoyancy and thus will 
have a feedback on the flow dynamics (as already discussed in Ref. 26)). It is be- 
cause of this interdependence of both instabilities that a dogmatic discussion about 
a dominant role of cither the one or the other at all times and for all conditions 
seems counterproductive and meaningless. 

A more serious possible argument against an important role of the SASI in real 
supernova cores at any time follows from the fact that all current 3D simulations that 
account for neutrino heating and cooling (Fig. 20; Refs. 18), 23), 24), 31), 34), 55), 70), 
73), 74)) do not show the development of large-amplitude dipolar shock oscillations 
as seen in 2D models. Even for the 27 Mq progenitor Ott et al.^^^ could not flnd any 
strong evidence of SASI activity once neutrino-driven convection has started (whose 
fast growth, according to the authors, is enabled by sizable numerical perturbations 
caused by the employed cartesian grid). 

Smaller saturation amplitudes of low-order shock oscillation modes in 3D have 
been explained by the distribution of turbulent kinetic energy over a larger number 
of available modes associated with the additional spatial degree of freedom in 3D 
compared to 2D.^^) This interpretation may be consistent with the observation that 
the direction of the dipolar deformation wanders stochastically™^ in 3D whereas its 
axis is artificially fixed to the polar symmetry axis in 2D models. Nevertheless, 
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the explanation is neither fuhy satisfactory nor generahy vahd, because it imphcitly 
assumes that the total kinetic energy stored in turbulent flows of the gas in the 
postshock layer, in particular also in flow structures on the largest possible scales, 
cannot become considerably higher in 3D than in 2D. A priori a justification of 
this assumption is not obvious, but for some reason higher kinetic energies of the 
postshock layer indeed seem to be disfavored as suggested by the results in Ref . 24) , 
which showed roughly the same nonradial kinetic energies in 2D and 3D models near 
the explosion. 

Although current 3D simulations are still handicapped by a large number of 
approximations in the employed microphysics and neutrino transport, and although 
they mostly study highly simplified, parametrized setups^^^'^^^'^^^'^'^) or suffer from 
very limited numerical resolution ^'^^ or artificial, grid-induced perturbations-'^^)'^*')''^^) 
(or even other, so far not realized, numerical problems in codes newly applied to the 
stellar core-collapse problem) , the basic similarities of the outcome of the increasing 
pool of simulations lends growing support to expectations that SASI is subdominant 
in 3D.^^) Instead, neutrino heating and associated buoyancy appear to be the main 
trigger of turbulent mass motions in the postshock layer and of the shock deformation 
seen in the numerical models. 

Reflecting current results of 2D and 3D simulations, the conclusion seems to 
be unavoidable that the turbulent dynamics of the postshock mass motions in 3D is 
different from 2D and that this is likely to also affect the way how the nonradial flows 
in the neutrino-heating layer exactly foster the transition of the accretion shock to 
runaway expansion. If fully self-consistent, more sophisticated, well resolved, and not 
artificially grid-perturbed 3D models will confirm the current results of simplified, 
incomplete, deficient, and parametrized 3D setups, in which various feedback effects 
are ignored, we will be ready to agree with Dolcncc et al."^^^ that the supernova 
mechanism in 3D operates differently from 2D. It may be expected that the basic 
features and implications of 3D turbulence as seen in current models will be retained 
even in more elaborate simulations (unless the numerical grid of the simulations has 
a nonneglibile influence). It is less clear, however, whether the quasi-stationary shock 
expansion, which is less time-variable and less dynamic than in the 2D models, will 
survive when all relevant feedback effects are taken into account. 

In any case, however, we do not think that these aspects imply a fundamental 
revision of the overall scenario of how core-collapse supernovae achieve to explode. 
In the end they only mean refinements in details, although, unquestionably, such 
refinements could be quantitatively relevant and even crucial. Nobody has rejected 
and can seriously deny the important role of neutrino-driven convection also in 2D. 
Independent of whether or not SASI has a supplementing function in 2D, whereas 
potentially being of reduced or minor relevance in 3D, the basic combination and 
interplay of ingredients in a working supernova mechanism will be same: Multi- 
dimensional hydrodynamic instabilities help pushing the shock to larger radii, thus 
stretching the dwell time of shock-accreted matter in the gain layer and enhancing the 
efficiency of neutrino-energy transfer. In the case of successful explosions, bubbles 
and large-scale plumes of neutrino-heated matter will rise (whether constrained by 
axisymmetry or not) to ultimately bring the shock up to the runaway threshold. The 
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mechanism is powered by neutrino-energy deposition, which also provides the energy 
input to the blast wave for unbinding the stellar envelope. The functioning of this 
neutrino-driven mechanism is crucially supported by hydrodynamic instabilities. 

Not each and every single dynamical aspect that needs revision or receives re- 
finement in the generalization from 2D to 3D modeling — as important for a precise 
and quantitative understanding it may be — will mean a fundamentally new twist 
of this basic picture of the mechanism. As an example in this context we mention 
the discussion of the runaway growth of buoyant, neutrino- heating-driven bubbles by 
Dolence et al.^^-* The elements of this picture are not generically linked to the 3D sit- 
uation (nor is this so for the discussion of the simplified, analytic toy model presented 
in the paper) , they also apply to 2D geometry, and the basic outcome of the analysis 
concerning shock evolution as function of time as well as critical luminosity condition 
remind one of similar results for the ID case (apart from a quantitative reduction of 
the threshold value of the critical luminosity by multi-dimensional effects). 

§6. Summary and conclusions 

We have discussed results of 2D simulations of the Garching group, in which 

for a growing set of progenitor stars successful, albeit mostly weak, neutrino-driven 
explosions could be obtained (Sects. 3 and 4). Most of these simulations were per- 
formed with a fully relativistic treatment and some with an approximate description 
of relativistic effects, but all of them included in a consistent manner the effects 
of hydrodynamics and of sophisticated, energy-dependent neutrino transport (cither 
with ray- by-ray approximation or two-dimensional, two- moment closure). Despite 
interesting and quantitatively relevant differences in many details, the successful runs 
with relativistic physics are basically compatible with the simulations of the same 
progenitors with relativistic approximations. 

We also attempted an assessment of these 2D results in the context of the current 
dispute about the dimensional dependence of the hydrodynamics of core-collapse 
supernovae (Sect. 5). Since 3D fluid dynamics in the supernova core seems to differ 
from 2D flows in a variety of aspects, what can one learn from self-consistently 
exploding 2D models? What can 2D modeling tell us at all, in particular since the 
explosions turn out to be weak and not to be able to reproduce observed supernova 
energies? 

First, the 2D simulations show that multi-dimensional effects in the most so- 
phisticated neutrino-hydrodynamics supernova models indeed lift the core conditions 
very close to the critical threshold for explosions, in most of the cases explored in 
Sect. 4 even beyond. This is obviously different from ID modeling, where failures are 
obtained for all but the lowest-mass supernova progenitors. We therefore consider 
the 2D successes as quickening progress and as an excellent reason to appreciate the 
glass of wine — measuring our understanding of the supernova mechanism — as half 
full. Moreover, the successful explosions in 2D are certainly a very important and 
necessary confirmation of the notion that the critical luminosity needed for runaway 
conditions is reduced in the multi-dimensional case. A confirmation by sophisticated 
and fully self-consistent models was unquestionably necessary, because all other an- 
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alytical and numerical studies that reached this conclusion were done with highly 
idealized setups, in which most of the complex feedback processes in real supernova- 
core environments, especially also dissipative effects, were ignored.^^'"^"^'"^^^'^^^'^^^ 

Second, one should keep in mind that the neutrino-driven mechanism is powered 
by neutrino heating, which therefore is the crucial ingredient. Hydrodynamic insta- 
bilities, as indispensable as they may be, play "only" a supportive role. Although the 
exact way they operate is, of course, important for our theoretical understanding and 
demands detailed explorations, dimension-dependent differences like the question 
whether SASI plus neutrino-driven convection or mere neutrino-driven convection 
and 3D turbulence are the assisting agencies of shock expansion, are possibly more 
a kind of refinement than an aspect of central nature in the fundamental problem of 
the explosion mechanism. Despite the constraint to two dimensions, successfully ex- 
ploding models therefore lend strong support to the viability of the neutrino-driven 
mechanism. 

The recent literature on the explosion mechanism in different dimensions con- 
tains a variety of incorrect claims. One of them, for example, is the argument that 
local neutrino heating associated with ray-by-ray neutrino transport unphysically 
correlates with low-order mode shock and matter motions and unphysically pumps 
them.^^) Generally, however, a problem of this kind appears highly unlikely because 
violent, bipolar motions of the shock and postshock layer were seen in 2D models 
with ray-by-ray transport as well as with true 2D transport schemes (cf. Fig. 17 
and, e.g., Refs. 53), 54), 69)) and even with simple neutrino-heating terms based on 
a spherical neutrino-lightbulb description (e.g., Refs. 14), 23), 24)). Moreover, vigor- 
ous sloshing motions of the shock were neither found for all investigated progenitors 
in 2D simulations (see Fig. 10) nor were they obtained in 3D models with ray-by- 
ray neutrino treatment (Fig. 20), although such shock dynamics was present in 2D. 
Ray-by-ray transport and large-amplitude SASI sloshing of the shock therefore do 
not seem to have any tight causal connection. Another incorrect claim concerns an 
alleged proposition by Hanke et al.^^^ that 3D simulations exhibit the tendency to 
become more similar to ID results. Hanke et al. reported their observations and 
partly speculative interpretation of results for a very specific, highly artificial and 
simplified modeling setup adopted from other works for comparison, in which they 
found less readiness for explosions in 3D simulations with increasing grid resolution. 
This result is still not understood in a broader context, but in the unsuccessful 3D 
cases the shock evolution and some, but not all, explosion-relevant quantities (natu- 
rally and undisputably) became more similar to ID simulations, and the postshock 
flow showed less vigorous mass motions than in the corresponding, successful 2D 
models. The statement by Hanke et al. was therefore made in reference to these 
special results but was not meant to be a prediction for 3D supernova dynamics in 
general. 

Current simulations in different dimensions suggest interesting possible differ- 
ences of the flow dynamics in supernova cores between ID, 2D, and 3D. In (artifically 
exploded) ID models the transition from shock stagnation and accretion to runaway 
shock expansion takes place either through a number of radial shock pulsations or 
by a continuous growth of the shock radius until accelerating expansion sets in.^^' 
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While self-consistent ID simulations do not produce explosions, 2D models show 
a reduced neutrino-luminosity threshold for the onset of explosions. In this case 
violent, nonradial shock oscillations (potentially connected to SASI activity) with 
increasing amplitude are seen for many but not for all progenitors. The conditions 
in the progenitor seem to play a crucial role for the type of observed behavior. More 
violent shock oscillations appear to occur in more massive stars with compact cores. 
Less massive progenitors with their less compact cores and faster decay of the mass 
accretion rate exhibit more signs of classical neutrino-driven buoyancy, in which case 
the shock expansion is more spherical and more continuous. In 3D, at least as far as 
current models of different sophistication — but only few with self-consistent coupling 
of hydrodynamics and neutrino transport,^^)'^^^'^^^'^^^'^^^ and all of them still with 
many simplications and deficiencies — allow one to conclude, the shock expansion 
appears to be more quasi-stationary, driven by the inflation of neutrino-heated bub- 
bles. The turbulent fragmentation of these bubbles leads to potentially large shock 
deformation but seems to push the shock outward in all directions more continuously 
than in 2D. 

While the first 3D results are inspiring and show new possible twists and direc- 
tions in the interplay of different physical ingredients of the highly complex problem 

of the supernova explosion mechanism, these first 3D simulations leave more open 
questions than they can answer. While general agreement exists about the supportive 
influence of multi-dimensional postshock flows in 2D, present results are inconclusive 
whether there is a further reduction of the critical luminosity in 3D compared to 2D. 
The results of Hankc ct al.^^^ arc in conflict with corresponding claims by Nord- 
haus et al.,^^-* and the more recent results of the Princeton group^^^' '''^^ also show a 
considerable weakening of the effect advocated previously: Although the newest 3D 
simulations of the Princeton group still yield somewhat earlier (~100-200ms) ex- 
plosions than their 2D models, the critical luminosity differences found before have 
nearly disappeared. A general trend to faster explosions in 3D, however, is neither 
clearly supported by the results of Ref. 24) nor by the outcome of more sophisticated 
models by L. Scheck.'''^^ So does 3D turbulence foster the onset of explosions even 
better than 2D flows? And if so, how large is the difference and by which effects is it 
caused? Turbulence in 3D pumps energy into smaller scales, but the energy stored 
on these scales is much lower than that on large scales. How can this have a consid- 
erable, supportive influence on the shock revival? Only a smaller fraction of the fluid 
elements (tracer particles) in 3D simulations^^-*' ^^-^ is found to possess longer dwell 
times in the gain layer, but a significantly larger fraction has shorter ones than in 
2D, which is fully compatible with the longer neutrino-heating timescales and lower 
net heating rates seen in 3D. If it is not more efficient neutrino heating, which 3D 
effects could account for faster explosions than in 2D models? Is it possible that 
different dimension-dependent effects in opposite directions compete with each other 
and partially compensate each other, which might explain seemingly contradictory 
trends of 2D vs. 3D found by different groups? Conclusive answers have not been 
given yet. In view of the turbulent energy cascade, are vigorous low-order modes of 
postshock mass motions generally excluded in 3D? Does neutrino-driven convection 
possibly destroy the coherence of SASI sloshing and spiral modes? Does SASI play 
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a role at all in supernova cores or are all relevant mass motions in 2D and 3D pure 
neutrino-driven convection as hypothesized in Ref. 35). Again, conclusive answers 
do not exist. If 3D turbulence on small scales is the crucial phenomenon that dis- 
criminates the onset of explosions in 3D from those in 2D, what is the influence of 
numerical resolution, what are the consequences of different grid treatments (spheri- 
cal vs. cartesian, adaptive meshes vs. fixed ones) and of different hydro solvers? And, 
on the level of physics, what is the influence of neutrino viscosity and in particular 
of magnetic viscosity? If small-scale turbulence is relevant, the amplification of the 
fields may not be negligible and may have an influence on the dissipation scale and 
on the transport of energy. 

In order to receive reliable answers, self-consistent models with all relevant in- 
gredients will be needed. Not every simplified approach that may be sufficient for 
studying observable consequences and measurable signals connected to supernova 
explosions initiated by the neutrino-driven mechanism, is equally justified when the 
explosion mechanism itself and its determining factors are to be explored. The lack 
of self-consistency and of important feedback effects may produce misleading results 
where subdominant effects show up much more strongly than they would if all rele- 
vant physics were included. In order to reliably account for the feedback of neutrino 
transport, generic transport effects have to be included and sufficient resolution is 
indispensable, in particular also in the surface-near neutron star layers where the 
density gradient steepens with time and the accretion-regulating cooling region is 
squeezed into a narrow shell with tremendously increasing cooling rates per nucleon 
(in Ref. 5) the cooling rates at late postbounce times were at least 10 times higher 
than those shown in Ref. 18)!). Also numerical artifacts potentially caused by grid 
effects and grid perturbations or numerical diffusion, e.g. connected to adaptive mesh 
refinement, will have to be investigated. Finally, if different hydrodynamic instabili- 
ties compete in their growth, the results might not be independent of the initial seed 
perturbations, and the still unsettled variations of physical quantities in the convec- 
tive regions of the progenitor cores might predetermine the growth of the instabilities 
after core bounce. ^^^'^'''^ Not only is a reliable definition of these initial conditions 
needed; also (noise-free) dynamical simulation codes that allow for a tracking of the 
growth behavior of these initial asymmetries during infall and postbounce evolution, 
will be necessary. 

Supernova modelers have only now begun to touch the vast wealth of dynamical 
phenomena that may play a role in the explosion mechanism of neutrino-heated 
supernova cores in the third dimension. Astonishing discoveries may be waiting for 
us when we continue to move forward into this unexplored territory. Don't let us 
freeze in humility in the face of this age-old challenge but let us get down to work! 
A lot of fun will be our reward! 
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